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ABSTRACT 



Cosmological Perturbation Theory (PT) is a useful tool to study the cumulants 
of the density and velocity fields in the large scale structure of the Universe. In papers 
I & II of this series we saw that the Spherical Collapse (SC) model provides with 
the exact solution to PT at tree-level and gives a good approximation to the loop 
corrections (next to leading orders), indicating negligible tidal effects. Here, we derive 
predictions for the (smoothed) cumulants of the velocity divergence field 8 = V • V for 
an irrotational fluid in the SC model. By comparing with the exact analytic results 
by Scoccimarro & Frieman (1996), it is shown that, at least for the unsmoothed case, 
the loop corrections to the cumulants of 8 arc dominated by tidal effects. However, 
most of the tidal contribution seems to cancel out when computing the hierarchical 
ratios, Tj = (8 J ) / (8 2 ) J ^ 1 . We also extend the work presented in Papers I & II to 
give predictions for the cumulants of the density and velocity divergence fields in non- 
flat spaces. In particular, we show the equivalence between the spherically symmetric 
solution to the equations of motion in the SC model (given in terms of the density) 
and that of the Lagrangian PT approach (in terms of the displacement field). It is 
shown that the ^-dependence is very weak for both cosmic fields even at one-loop (a 
10% effect at most), except for the overall factor /(fi) that couples to the velocity 
divergence. 

Key words: cosmology:large-scale structure of Universe-cosmology: theory-galaxies: 
clustering- methods:analytical. 



1 INTRODUCTION 

One of the aims of modern Cosmology is to understand the 
origin and growth of the large scale structure as we observe 
it today. Cosmological Perturbation Theory (PT) has re- 
vealed to be a key analytic tool to make predictions about 
the cosmological matter density 5 and velocity V fields in 
the weakly non-linear regime (i.e., when v,8 < 1). The first 
analytic results obtained were only applicable to the un- 
smoothed fields, Gaussian initial conditions (GIC) and flat 
space (i.e.,0, = 1 for a vanishing cosmological constant, 
A). These results concern the tree-level (leading-order for 
GIC, see Peebles 1980, Fry 1984, Bernardeau 1992, here- 
after B92). Recently, loop corrections to the tree-level am- 
plitudes (higher-order contributions for GIC) were obtained 
within the diagrammatic approach by Scoccimarro & Frie- 
man (1996a, 1996b) although their results basically concern 
the unsmoothed fields. Comparison with N-body simulations 



and observations made it necessary to incorporate the ef- 
fects of smoothing in the PT calculations (Goroff et al, 
1986, Juszkiewicz et al., 1993, Bernardeau 1994a, 1994b, 
Baugh, Gaztanaga, Efstathiou 1995, Gaztanaga & Baugh 
1995). Further results were obtained for non-Gaussian IC in 
a number of papers (Fry & Scherrer 1994, Chodorowsky & 
Bouchet 1996) which allows to discriminate among models 
of structure formation (Silk & Juszkiewicz 1991, Gaztanaga 
& Mahonen 1996). 

In Paper I (Fosalba & Gaztanaga 1998) and Paper II 
(Gaztanaga & Fosalba 1998) of this series it was shown that 
the Spherical Collapse (SC) model gives the exact tree-level 
contribution to PT irrespective of the nature of the IC and 
the geometry of the universe. This means that, when it comes 
to computing the one-point cumulants at tree-level, only the 
spherically symmetric solution to the equations of motion 
is relevant (see Paper I §3.6). This solution generates the 
monopole term, v n , i.e., the angle average of the kernels in 
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Fourier space. As a result, the evolved density field is entirely 
determined by a local transformation of the linear density 
field alone, what we shall call a local- density transformation 
(see Paper I §3.4), 



s = M) 



nl 

71 = 1 



['V 



(1) 



Making use of this property we were able to present in Pa- 
pers I & II results for Gaussian (GIC) and non-Gaussian 
IC (NGIC) respectively, for both the unsmoothed and 
smoothed fields (for a top-hat window). Furthermore, the 
SC model was used to make predictions beyond the tree- 
level in both cases which were shown to be in good agree- 
ment with N-body simulations up to the scales where PT 
itself is expected to break down (see Paper I §5 and Paper 
II §3.3). 

In this paper, we derive results for the peculiar velocity 
field V which is the counterpart to the density field in the 
Newtonian equations of motion. The peculiar velocity field 
as traced by a population of galaxies has the potential ad- 
vantage, with respect to the corresponding matter density, of 
being a better (or at least different) tracer of the correspond- 
ing underlying field. This is because galaxy velocities should 
follow the gravitational field regardless of their luminosity. 
These property has been proposed as a way to make pre- 
dictions about the density parameter of the universe Q (see 
Bernardeau 1994b) and references therein). In what follows 
we shall assume that the fluid is irrotational as is usually 
done in PT. This is a good approximation in the SC model 
as the nonlinear density field is given in terms of the lin- 
ear field alone, for which, PT predicts a dilution of vorticity 
with the expansion. Within this approximation the velocity 
field is entirely determined by its divergence, 8 = V ■ V. 
The cumulants (or connected moments) of the velocity di- 
vergence are defined as, {8 J ) C and the corresponding hier- 
archical ratios are Tj = {8 J ) c /(9 2 ) J_1 . Results for the one- 
point cumulants of the velocity divergence field were given in 
Juszkiewicz et ai, (1993) for the skewness T3 (at tree- level 
and GIC, alone). Bernardeau (1994a, 1994b, B94a,b here- 
after) systematically derived the dependence of these statis- 
tical quantities on the initial power spectrum for a top-hat 
window function and the strong (overall) omega dependence 
(see B94b) to leading-order (tree- level) for GIC. 

We shall also extend the work to non-flat spaces 
(i.e.,0, 7^ 1 for a vanishing cosmological constant A) for 
the cumulants of both the density and velocity divergence 
fields. Curvature effects on the hierarchical amplitudes have 
been investigated in a number of papers: Martel & Freudling 
(1991), and Bouchet et ai, (1992) obtained the (weak) Q 
dependence of the skewness of the density field in the La- 
grangian PT approach. Bernardeau (1994b, B94b hereafter) 
numerically integrated the ^-dependence of the skewness 
and kurtosis of the density and velocity fields. Exact ana- 
lytic results , concerning the curvature-dependence of the 
skewness, were derived by Bouchet et al, 1995 (BCHJ95 
hereafter) in Lagrangian PT and Catelan et al, 1995, in 
Euler space. Recently, Nusser & Colberg (1998) have inves- 
tigated the dependence on the cosmological parameters Q 
and A in the equations of motion and give some results for 
the variance and the skewness. 

This paper is organized as follows. Section §^ presents 



results for the velocity divergence field while in §|3j results for 
the ^-dependence in the SC model are given. A discussion 
and the final conclusions are given in §^j. Appendix §|X| pro- 
vides the connection between the vertex generating function 
in PT and the SC equation for arbitrary Q,. Appendix §|§| 
gives an alternative derivation of the Q dependence in the 
monopole approximation to Lagrangian PT. 



2 RESULTS FOR THE VELOCITY 
DIVERGENCE FIELD 

In this section we apply our results from the SC model for 
the density field (see Papers I & II) to derive the one-point 
cumulants for the velocity field in the quasi-linear regime. 
For the current analysis, we shall consider an irrotational 
fluid and, as a result, the velocity field is fully described by 
its divergence, 8 = W/H, where TL(t) = a(t)H(t) is the 
conformal Hubble parameter and a(t) is the scale factor. 
As mentioned above, this is a good approximation in the 
SC model as nonlinear dynamics are fully determined by 
the linear growth of perturbations for which vorticity decays 
with time. 

The velocity divergence is locally related to the linear 
density field through the continuity equation, which for an 
Einstein-deSitter space reads, 



8 + (1 + 5)8 =0, 



(2) 



where the dot denotes total time derivative. It follows that 
the linear solution is 81 — —Si. This way, introducing Eq.Jl] 
above, we can obtain an analog local-transformation for the 
velocity divergence, say, 



(3) 



where the unsmoothed fj,k coefficients are to be derived or- 
der by order in a perturbative expansion of the continuity 
equation. Thus one finds, 
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and so on. 



2.1 Smoothing Effects 

The smoothing effects for the velocity field for a top-hat 
window can be easily derived by applying similar argu- 
ments to those used for the density field (see Paper I §4.4 
and B94b). In particular, the continuity equation ensures 
that the smoothed velocity divergence field is proportional 
to the unsmoothed one, 8 ~ 9 as induced by the sharp 
cut in the smoothed density field with respect to the un- 
smoothed one, 8 ~ 8. Moreover, the local character of the 
unsmoothed velocity divergence field (see above) tells us 
that 8 ~ g[8i]- Now, from the linear continuity equation 
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we get, §i = —Si = -{ai/ai) Si = (<ri/&i) 9i. Introducing this 
last equality in the local transformation g which determines 
9 we finally get, 



(5) 



where ~ just means that the transformation is performed in 
Lagrangian space so that it has to be properly normalized 
when going back to Euler space (see Paper I §4.2). 

For a power-law power spectrum P(k) ~ k n , one finds 
(see Paper I), 



oi = ai (1 + 6) 



-7/6 



(6) 



where 7 = -(n + 3). As a result, the smoothed non-linear 
velocity divergence field is given by, 

§[§ l ]^e[§ l (i + s)' l/ % (7) 

to be normalized when transforming back to Euler space 
in an analogous way to the density field (see Paper I, Sec- 
tion 4.4). Notice that these results are formally equivalent 
to those provided by B94a for the smoothed fields in Euler 
space although, here, we work in Lagrangian space. 

According to this transformation, the first smoothed 
coefficients of the velocity field read, 



= 



Ma 



1 

4 



/'4 
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27 
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(8) 



and so forth, where in the above equalities we have applied 
the identities given by Ecl^. 

Combining Eqs.||,[^|] and we can derive the 

smoothed cumulants for the velocity divergence up to an ar- 
bitrary perturbative order in an Einstein-deSitter universe 
for a top-hat window and a power-law power spectrum, for 
QIC & NGIC. 



2.2 Gaussian Initial Conditions 

Following the notation introduced in Paper I §5, we shall 
use for the perturbative expansion of the variance, 



2 



= (§*) = 



, s 2 ,i °i 



(9) 



where S24 = 1 and the subscript in the coefficients labels the 
order of the perturbative expansion. Note that a; denotes 
the rms fluctuation of the linear density field. As for the 
hierarchical amplitudes, we keep the above notation with the 
added labeling of the order of the moment, J, that defines 
the Tj coefficients, 



Tj 



<* 2 > J ~ 



Tj,i a 1 . 



(10) 



As we did with the density field (see Paper I §5), we shall 
denote the leading order contributions (i.e., the tree- level for 
GIC) by Tj 0) . 
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Table 1. Values for the higher-order perturbative contributions 
in the SC model for the unsmoothed (n = —3) and smoothed 
(n = —2,-1,0) velocity fields. 



2.2.1 Spherical Collapse Results 

For GIC the odd terms in the perturbative expansion vanish, 
thus, 
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2,6 4 . B 6 , 
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(11) 

For the variance, skewness and kurtosis of 9 [for a top-hat 
window and a scale- free power spectrum, where n = — (7 + 
3)] in an Einstein-deSitter universe, we find: 
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1134 ' 648 
The results for Tj,o above are identical to those systemati- 
cally obtained by Bernardeau (1994a, 1994b) who made use 
of the vertex generating function formalism, although here 
they are restricted to a power-law power spectrum. They can 
be straightforwardly extended to an arbitrary (more realis- 
tic) power-spectrum by applying analog expressions to those 
for the density field provided in Appendix Al of Paper I. 
They also extend those results to the one-loop contribution 
within the SC approximation. Table |l| displays the results 
for the smoothed field for different values of the spectral 
index. 

To compare our results to those available in the liter- 
ature we should focus on the unsmoothed one-point cumu- 
lants, 

0$ fit of + 10.97 af + 0(af) 
-T 3 w 3.71 + 3.86 af + 0(af) 
T 4 w 27.41 + 131.06 a-, 2 + O(of). (13) 

The results from the exact PT in the diagrammatic approach 
derived by Scoccimarro & Frieman (1996a) are, 
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SC: 2nd. Order 
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Figure 1. Departures from the tree- level up to one loop as a 
function of the linear mis fluctuation. The variance, skewness and 
kurtosis of the velocity divergence predicted by the SC model 
are shown for different values of the spectral index: the dotted 
line shows the n = —3 (unsmoothed) case, and the short-dashed 
(n = —2), long-dashed (n = —1), dot short-dashed (n = 0), dot 
long-dashed (n = 1) depict the behavior for the smoothed field. 
The solid line shows the tree-level values. 



a 2 g ^af + 1.08 af + 0(af) 



(14) 



and the first corrective term for T3 was viewed to strongly 
depend on the spectral index due to non-local effects: 
-2.58 > T 3 , 2 > -4.64, for 2 > n > -2 . Thus, our lo- 
cal value T3,2 = —3.86, is compatible with the average of 
the above non-local values. For the variance however, the 
corrective term seems to be dominated by non-local (tidal) 
effects which cut down non-linearities in an order of magni- 
tude, unlike the case of the density field (see Paper I §5.1). 
The situation might be different for the smoothed cumu- 
lants, as we have seen (in Papers I&II) that for n ~ —1.5 
tidal effects seem to cancel out. 

Fig [l] depicts the departures from the tree-level contri- 
butions for different values of the spectral index as the linear 
variance grows. Notice that, in general, the smoothing effects 
tend to suppress nonlinear contributions in accordance with 
the trend observed for the density field (see Paper I Table 2 
& Paper II Table 1). 



2.2.2 Spherically Symmetric Zel'dovich Approximation 
Results 

The Zel'dovich Approximation (ZA hereafter, see Zel'dovich 
1970) formally corresponds to the first order in the La- 
grangian approach to PT (see Appendix |b| for details) which 
is based in a perturbative expansion of the Jacobian that re- 
lates the Lagrangian to the Eulerian coordinates. It assumes 
that the particle positions in comoving coordinates follow 
straight trajectories in Lagrangian space. In this section we 
proceed and give the results for the spherically symmetric 
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Table 2. Values for the higher-order perturbative contributions 
in the SSZA for the unsmoothed (n = —3) and smoothed (n = 
-2, -1, 0) velocity fields. 



solution to the equations of motion within the Zel'dovich 
Approximation (SSZA hereafter, see Paper I Appendix A4). 
In particular, we give estimates of the cumulants for the 
smoothed velocity divergence (for a top-hat window and a 
power-law power spectrum) in an Einstein-deSitter universe, 
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thus, the unsmoothed cumulants (7 = 0) show the following 
scaling, 



o-q w at + 6.56 at + 0(af) 
-T 3 w 2 + 2.07 af + 0(<jf) 
Ta ~ 8 + 38.44 af + 0(af). 



(16) 



This should be compared to the unsmoothed results from the 
exact PT in the diagrammatic approach up to the one-loop 
contribution (see Scoccimarro & Frieman 1996a), 

aj^af + 0.73 af + O(af) 
-T 3 w 2 + 1.64a ; 2 + O(af) 

T 4 w 8 + 17.58 a 2 + O(af) (17) 



illustrating again the point made above (see §2.2.1): tidal 
effects dominate the variance of the velocity field and make 
non-linear effects small, what is at variance with predictions 
from the local SSZA dynamics. Again here, the situation 
might be different for the smoothed cumulants, as we have 
seen (in Papers I&II) that for n ~ —1.5 tidal effects seem 
to cancel out. Despite this limitation intrinsic to our local 
approximation, we see that tidal effects almost cancel their 
contributions to the reduced cumulants. This cancellation 
renders the SSZA approximation a good estimator of these 
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Figure 2. Same as Fig [l] for the SSZA. Note that the line n = 
— 1 is not shown for T3 and T4 since the tree-level hierarchical 
amplitudes vanish. The values n = —3 (unsmoothed) and n = — 1 
lead to the same scaling properties for the skewness and kurtosis 
so the corresponding lines overlap. 



quantities within the ZA dynamics, in line with the results 
obtained previously for the density field. Of course the re- 
sults in EqJlS] are a better approximation to the exact PT. 

Table fe shows the results for the smoothed velocity di- 
vergence field in the SSZA for different values of spectral 
index. Fig ^| displays deviations from the tree-level contri- 
butions for different values of n as the linear variance ap- 
proaches unity. 



2.3 Non-Gaussian Initial Conditions 

For NGIC, no additional equations need to be solved as all 
the information relevant for the computation of the one- 
point cumulants is encoded in the Gaussian tree-level (see 
Paper I §3.8). The perturbation expansion is worked out in 
an analog way to that for GIC, the only difference being 
the order at which the dominant terms appear, as induced 
by the scaling properties of the IC (see Paper II §4). We 
represent the different orders in the expansion as we did for 
the GIC above, Eq.§-@. 

We will assume, as a generic case, that the NGIC follow 
a dimensional scaling: {0f) c oc {0j) J ^ 2 , and we will define the 
corresponding initial amplitudes as: 



Bj 
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not to be confused here with the corresponding coefficients 
for the density field. This scaling induces different ordering 
in the perturbative series as compared to the GIC one. In 
particular, we find the following results for the first contri- 
butions to the cumulants filtered with a top-hat (assuming a 
power-law spectrum) and an Einstein-deSitter background, 
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where T G correspond to the tree-level PT results (7j,o) for 
GIC. The unsmoothed one-point cumulants (setting 7 = 
in the Gaussian coefficients above) we find: 
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We can compare these results to those from exact PT, de- 
rived by Protogeros and Scherrer (1997), which give for T3 
for the unsmoothed fields (setting the overall S7-dependence, 
/(fi) = 1 for a flat universe), 

10 
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being J[^j] = D(t) J I[£j], where D(t) is the linear growth 
factor and the I[£j] denote some integral of the initial J- 
point functions which are intrinsically non-local terms (see 
Eq.[30] in Paper II). Thus, using T 3 G = -26/7 and Bj = 
{9/)c/o-f in our expressions (see Eq.Ji^] above), we find, as 
we did with the density field (see Papers I & II), that the 
SC model gives the exact result up to some tidal (non-local) 
terms which result in a small contribution to the hierarchical 
amplitudes. This smallness is partially due to the fact that 
they enter as a difference between integrals of the initial J- 
point functions which are generically of the same order, at 
least for the density field (see Paper II, Section 4.4). 

For the smoothed fields however, our results must be 
taken as a prediction since there are no explicit results for 
non-Gaussian initial conditions. For some general implicit 
results, which depend on the initial J-point functions for the 
smoothed T3, see Protogeros & Scherrer (1997). This should 
contain all our local terms when given in a somewhat more 
explicit fashion. 

Figs |^ and ^ show the evolved variance, skewness and 
kurtosis as they depart from their first perturbative contri- 
butions rj 0) . As we pointed out in the case of the density 
field (see Paper II, Section 3.1, Figs. 1 & 2), the variance of 
the velocity divergence is seen to be dominated by the ini- 
tial conditions in the second perturbative contribution (only 
terms ~ Bj appear), while gravity takes over once the third 
perturbative contribution is included. This is realized by the 
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Figure 3. Same as Fig |l| for non-Gaussian dimensional initial 
conditions. We assume Bj rj 1. 



Figure 4. Same as Fig ^ when the 3rd. perturbative contributions 
are taken into account in the evolved one-point cumulants. 



enhancement of the non-linear variance with respect of the 
linear one (at least, for large values of the linear rms den- 
sity fluctuation) in accordance with the expectation from 
the SC model for GIC. The hierarchical amplitudes show 
the same coupling to the initial conditions as for the density 
field so all perturbative contributions are taken over by the 
initial conditions. In Figures ^| & ^ we have taken the ini- 
tial dimensional coefficients B.j ~ 1, in analogy to the case 
for the density field (see Paper II §3.1), although we find 
no physical motivation for why they should be comparable. 
However, we have checked that the qualitative behavior of 
the statistical quantities and thus the overall picture is not 
significantly altered if we change those amplitudes by an 
order of magnitude i.e., our conclusions are robust. 



3 THE ^-DEPENDENCE IN THE SPHERICAL 
COLLAPSE MODEL 

3.1 The fi-dependence in the Density Field 

3.1.1 The Spherical Collapse Model Approach 

In the previous papers (see Papers I & II) we assumed a 
spatially flat (Einstein-deSitter) universe with a vanishing 
cosmological constant A, in order to derive our predictions 
for the one-point cumulants. Although it is the simplest an- 
alytic approach and there are strong theoretical arguments 
which favor an Q. = 1, A = universe {e.g., standard infla- 
tionary models), it remains to be seen whether or not this 
is actually so. 

In our calculations, we have neglected the effect of a 
non-vanishing cosmological constant A. However, detailed 
analyses have shown that the dependence of the cumulants 
on this parameter is extremely weak anyway, at least for a 
flat universe O + A = 1 (see Lahav et al, 1991, BCHJ95). 

We showed in Paper I that in the general case of a 



perturbative approach to estimate the leading order (tree- 
level) contribution to the cumulants, it is not necessary to 
calculate the full kernels F n , but only the monopole contri- 
bution, v n =< F n >, which corresponds to the spherically 
symmetric (angle) average. Thus the values of v n can be 
determined by finding the spherically symmetric solution to 
the equations of motion. We also show that for the gravi- 
tational evolution of the density field, this solution is given 
by the spherical collapse (SC). The proof (§4.1 in Paper I) 
holds irrespective of the geometry of the universe, i.e., for 
arbitrary f2 (which shows explicitly through the Friedmann 
equation: AnGp = 3/2QH 2 ). Therefore, we can write the lo- 
cal non-linear transformation introduced in flat space (see 
Paper I, §4 Eq.[36]) and write, 



E 



[Si. 



(22) 



where the v n (Q.) coefficients are a generalization of the flat- 
space monopole amplitudes, u n , to non-flat FRW geome- 
tries. This can be readily obtained by considering the para- 
metric solution to the SC dynamics in an open universe 
($1 < 1): 



DM 



sinh 6 - 6 



sinh ip — ip 



2/3 



/ cosh ip — 1 V 
V cosh 9 - 1 ) 

for Si > —D(ip), and 



sinh 9-6 
sinh ij) — ip 



5i = -D(i>) 



- sin 9 + 9 
sinh ip — ip 



2/3 



+ 1 



cos ib — 1 \ 3 / sin 9 — 9 
1 x ' 



- cos 9 + 1 



\ sinh ip — ib 



(23) 
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for 5i < —D(ip), with 



D(1>) 



9 sinh?/>(smhi/! — tp) 
2 (cosh Tp - l) 2 



= cosh 



(24) 



where D(%j)) is the growing mode of the linear density con- 
trast (see also B92). In order to get the solution for an Q > 1 
universe, the transformation ip — > iip must be performed 
in the above expressions. This way, the v n (Sl) coefficients 
may be numerically integrated by expanding the paramet- 
ric solutions around Si = 0. An alternative derivation of the 
f2-dependence of the Gaussian tree-level was presented by 
Bernardeau (1992) in terms of the vertex generating func- 
tion. In Appendix Al we present the connection between the 
latter and the SC equation for arbitrary Q. 

A good fit to the above results can be obtained by re- 
alizing that in Lagrangian PT (see §|b|) the ^-dependence 
of the hierarchical amplitudes only appears beyond the first 
perturbative contribution (which corresponds to the ZA). 
As the larger the density parameter is, the smaller the Sj 
parameters are (this is a monotonic behavior, at least within 
the range 10 > il > 0), one is led to the simplest possible 
Ansatz for Sj(Q), which is of the form, 

Sj(Q) = Sj A + BQ- a 

being a and B positive constants. Imposing that for Q — 1 
one must recover the well-known Einstein-deSitter values, 
we find: 



Sj(n) = s z j A + (s: 



j 



s z j A ) n- 



(25) 



so that a is the only parameter to be fitted by inspection 
in a log-log plot of the numerically integrated solutions. In 
particular, for the cosmologically favored range 2 > Q. > 
0.1, we obtain the following numerical fits for the tree-level 
hierarchical amplitudes, 



S 3 (fi) 
S 4 (f2) 



4i/ 3 (fi) + 12 v 2 W 2 
60712 20728 



■Q 



-3/91 



1) 



- 1323 ' 1323 i° ' ^ ™ 

from which we derive the fits for i^(fl) and z/3 (Q). Replacing 
these unsmoothed coefficients in the general expressions for 
the smoothed Vk (since the smoothing effects were derived 
regardless of the geometry of the universe, see Paper I §4.4) , 
we get at tree level, 

Ss(fi) = 3^(?7)^y + ^ (fT 3/91 - 1) + 7 

-2 



S 4 (fi) = 4*s(fi) + 12 1/2 (n) 



60712 20728 



1323 



+ 



1323 



'n- 



-1/27 



-1 



, 7 2 , , / -3/91 
y7 + g7 +47 (^ 



(27) 



Moreover, the numerical fits for the v n (Sl) can be used 
to derive the corresponding approximate ^-dependences for 
the cr-corrections to the cumulants as they carry all the dy- 
namical information in the SC model (see Paper I, §3.5). 
Fig |i] shows the departures from the tree-level amplitudes 
for different values of Q. As shown there, the dependence is 



rather weak and only for values of Q close to zero are the 
predictions a 10% greater than those for flat space. 

We have also explored the f2-dependence of the non- 
linear variance. Making use of the fits derived for the i/fc(fi) 
above, we find for the first non-linear correction (one-loop), 



S2,4 



1£2£+143 +11 1° to" 6 / 91 - 
1323 126 ' 36 ' 49 v 

fn- 1 / 27 _ 1) + JL (n 7 - 64) («- 3 / 91 



1323 



1) .(28) 



From the above result we estimate that the non-flat 
(fl < 1) variance to one-loop can be a 10% greater at most 
(for <7; < 1), than the flat space value, similarly to what we 
found for the Sj amplitudes. We stress nevertheless that the 
SC result for the variance only gives an accurate approxi- 
mation to the exact PT result (and the N-body predictions) 
around n ~ —1.5 (7 ~ —1-5), where tidal forces yield a 
negligible contribution (see Paper I §5.2). 

The general trend observed in the hierarchical ampli- 
tudes for the density contrast in the non-linear regime, ac- 
cording to the SC model, is that, the lower the density pa- 
rameter is, the higher the amplitudes are. This reflects the 
fact that non-linear terms in the perturbative expansion of 
the overdensity in the low SI model are larger at a given 
time than the corresponding ones in flat space. This seems 
in qualitative agreement with the results recently found by 
Nusser & Colberg (1998) concerning the variance and the 
skewness from N-body simulations, for standard CDM mod- 
els (see their Table 1). A direct comparison is difficult to do 
in a significant way given that there are no errors, which are 
important for this slight deviations, and the quoted scales 
lie close to the nonlinear regime a > 1. 



3.1.2 The Lagrangian PT Approach 

An independent derivation of the approximated in- 
dependence in the hierarchical amplitudes of the density and 
velocity fields is provided by the Lagrangian PT approach 
(see Moutarde et al., 1991, BCHJ95). Since the spherically 
symmetric solution to the equations of motion gives the ex- 
act tree-level contribution to PT (see Paper I §3.4 & Paper 
II §2.4), we shall impose the spherically symmetric condition 
on the relevant equations of motion in Lagrangian PT to get 
the exact tree-level as well in that formalism. Lagrangian PT 
is formulated in terms of the perturbation of the original 
trajectory in Lagrangian space q, so that the perturbative 
series is built by expanding the Jacobian J = \dx/dq\ of the 
transformation between Lagrangian and Euler coordinates, 
x(t,q) = q + ^(t,g) (being the so-called displacement 
field). This is equivalent to expanding the density contrast 
in Lagrangian space since, 

5=P--1 = \-1. 



J 



(29) 



In this framework, the derivatives of the displacement 
field, "Jij, play the role of the smallness parameter in the 
perturbative expansion, instead of the linear density con- 
trast, Si, as it was the case for the SC model. By imposing 
spherical symmetry in the growth of perturbations we set 
^1,1 = Vl/2,2 = ^3,3 = ^k,k (being zero otherwise). So for 
all perturbative orders, it holds that, J (n) tx [J (1) ] n (see §|§| 
for details). This is the key point in the connection with the 
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m 



in 

CO 

m 




a, 



Figure 5. Predictions for the departures from the tree- level con- 
tributions to the skewness and kurtosis for non-flat universes. The 
solid line corresponds to the fiat-space prediction while the short- 
dashed, long-dashed and dot-dashed lines give the behavior for 
models with Q = 0,0.4,2 respectively. It is seen how the weak 
dependence on with the density parameter of the non-linearities, 
makes it hard to distinguish them from those of flat-space un- 
less Q approaches zero. The upper group of lines corresponds to 
n = —3, while lower groups of lines show the smoothed cases: 
from n=-2 to n=l (bottom group of lines). 



SC model in the perturbative regime as it explicitly shows 
that all the perturbative orders can be expressed as powers 
of the linear one as was the case for the expansion in terms 
of the linear density field in the SC model. In other words, 
we can build a local-density transformation (see Eq.Jl]) in 
Lagrangian PT from the spherically symmetric solution to 
the equations of motion whose coefficients (the Qk functions, 
see §[§]) carry the same information than the Vk ones in the 
SC model. Identifying the linear terms in both expansions 
we can formally relate all higher perturbative orders one by 
one, from which we can derive the relations between the co- 
efficients (gk with i/fe) that determine the one-point statistics 
in both formalisms. 

In particular (See Appendix A2), we get for the first 
coefficients, 



- - § »-3 



9 V g\ gl gl 



(30) 



Once we relate both unsmoothed coefficients Vk with g k , 
we can simply apply all previously given results for the SC 
model by rewriting the expressions appropriately in terms of 
the new gt coefficients. It is important to stress the equiva- 
lence between the SC model predictions in Lagrangian and 
Euler space at tree-level (see Paper I §4.2). Thus only be- 
yond tree-level (next-to-leading terms in PT) there are dif- 
ferences between the SC model formulated in Euler space 
and the spherically symmetric approach to Lagrangian PT. 



Furthermore, the equivalence between both formalisms 
(the SC model in the perturbative regime and the La- 
grangian PT) in Lagrangian space, allows to induce a fit 
to the fi-dependence of the gt functions from that of the 
Vh ones. We can compare the fits given by BCHJ95 to our 
fits for the cumulants to see which ones match better the 
numerical results. Although the exact analytic expressions 
for gi & <?2 are known (see BCHJ95), no exact results are 
available for higher- orders, so we shall focus on their numer- 
ical fits for the cosmologically favored parameter space for 
f2 and work out the predictions. According to BCHJ95, near 
= 1, good fits are provided by the following factors, 



i-l> 



-2/63 gSa 
' 3 

9i 



-4/77 



536 

gf 



21 



-2/35 



gl 7 gf 3 

Combining Eqs.]3c|] & |n| we get the fits to the co- 
efficients from the Lagrangian PT approach. The smooth- 
ing effects for a top-hat window can be incorporated as 
described in Paper I, Section 4.4. The resulting smoothed 
coefficients Uk(Q) in terms of the unsmoothed ones i/fc(f2), 
are given by Eq. [A3] in Paper I. Replacing them in the 
smoothed Sj amplitudes yields, for a power-law spectrum 
> + 3)], 



(31) 



[with 7 = 

S 3 («) = 3W 
S 4 (fi) = 4T^+12T^ 



34 

y 



+ 



in 



-2/63 



60712 
1323 

(O 



+ 



48 
49 

4/77 



-4/63 



1) +7 

, 80 
> + 63 



-2/35 



+ 



27 
62 
3 



7 



7+3-7" 



(f + (32) 



Fig [] shows the ^-dependence of the skewness and kur- 
tosis of the density field at tree-level, which is shown to be 
very weak, as expected. In particular we compare the exact 
solution for the f2-dependence to the fits to the exact behav- 
ior according to the SC model approach (based on Eq.[p5t) 
on one hand, and those based on the spherically symmetric 
solution to the Lagrangian PT (for which we make use of 
the BCHJ95 fits to the g k functions, Eq.Q). Our fits are 
slightly more accurate than those given by BCHJ95, which 
nevertheless only underestimate the exact ^-dependence of 
the skewness, 53, and kurtosis, S4, by 3% and 5% at most. 



3.2 The f2-dependence in the Velocity Divergence 
Field 

In this section, we extend the results obtained in §2 to the 
case where the density parameter 57 7^ f. In general, an over- 
all factor /(fi) = d log D(t)/d log a(t) (see Eq.^j below), 
couples to the velocity divergence in the continuity equa- 
tion for a non-flat space -what we shall call the strong £7 
dependence-, which scales in the moments of the velocity 
divergence (9 J ) oc /(f2) J . Therefore, for the cumulants, one 
expects 



Tj(0) 



■Tj(n = i), 



(33) 



with /(f2) w fi ' , provided one assumes a vanishing cosmo- 
logical constant A = (see Peebles 1980). Nonetheless, there 
is also a weak dependence on the density parameter that 
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Figure 6. Comparison between the exact (numerically inte- 
grated) Sj(Q) amplitudes (normalized to the flat space values) at 
tree-level (solid line) and the numerical fits to the f!-dependence 
for open models (dashed lines). The case n = — 1 is displayed. 
The short-dashed line shows the SC predictions making use of 
the BCHJ95 fits to the functions, Eq. Jjl|], in the context of 
the Lagrangian PT, while the long-dashed line gives the fits to 
the SC model solutions according to the Ansatz Eq. |25| . 



modulates the dominant (strong) dependence (see BCHJ95 
and B94a). This dependence is induced by the density field 
in a non-flat space and is obtained through the equation 
of motion when the generalized coefficients of the density 
field Vk{Q) replace those for flat-space in Eq. [^|. However, as 
commented above (see §3), the fits to the z/fc(Q) coefficients 
can be derived either in the SC model approach or in the 
Lagrangian PT approach since they are formally equivalent 
(see Appendix ^ for details) . 

From the fit s to the SC model approach for the density 
field (see § ^.1 we induce the following approximate in- 
dependences for the velocity divergence field whenever 2 > 
Q > 0.1, 



-T 3 /(ft) = 3R^)- f + y (^ 3/91 



!)+■ 



T A f(ny = 4 M3 (^) + i2jia(n) 

w 4 



3022 + 12 (fi _ 6/91 



13 
7 



441 49 

7 (n~ s/91 - 



i 



, 169 7 2 
1 ) + l2- 7+ 12 7 

(34) 



Alternatively, when we make use of the VkiP) as they are 
derived from the Lagrangian PT approach (see Appendix |b|) 
we get, 



Ti /(fi) 2 



25 + H (n-V63 

7 7 v 



-1 +7 



12088 192 



441 

I (n- 



+ 



49 



4/77 



-4/63 



338 



-1 



80 / -2/35 

21 \" 



2T 7+ 3 7 



Figure 7. Same as Fig |ij for the weak S7-dependence of the ve- 
locity divergence (i.e., excluding the overall ~ /(f2) factor). 



+ 



- (22 + 13 7 )- 



in 



-2/63 



-1) 



(35) 



where the strong ^-dependence, /(ft) = d\ogD(tp)/d\oga 
(where a is the scale factor and D(ip) is the growing mode of 
the linear density contrast), can be exactly integrated and 
yields for fl < 1, 



3 sech(tp/2) ((2 + Q) tp - 6 VT^) 
Q (9 sinh(V>/2) + sinh(3^/2) -6V cosh(^/2)) 
while for Q > 1, it has the form, 



m) = 



3 sec(ifi/2) ((2 + 0) ip - 6 y/tT- 



■(36) 



n (9 8111(^/2) +sin(3V>/2) - 6ip cos(i/>/2)) 

where, tp = cosh _1 (2/fi - 1) . However, the fit f(Q) » Q, 3/5 
is a good approximation for open models as originally put 
forward by Peebles (1980). Other more accurate fits are 
/(ft) w fi 13/22 for open models and /(ft) w fi 13/23 for closed 
models within the range 2 > f2 > 1. 

In Fig 7 we plot the weak f2 dependence (i.e., excluding 
the overall ~ /(ft) factor) for the hierarchical amplitudes 
at tree- level, for a value of the spectral index n = — 1. It 
shows two independent fits, the Lagrangian PT one (based 
on the spherically symmetric approximation to the pertur- 
bative expansion of the displacement field) and those within 
the SC model approach (same approximation but applied 
to the expansion of the density contrast), against the exact 
behavior numerically integrated. It shows that the BCHJ95 
fits to the <7fc functions (see Eq.|3l)) systematically underes- 
timates the exact (weak) ^-dependence of the hierarchical 
amplitudes at tree-level, T3/ and T4/ 2 , for open models (in 
the range 1 > ft > 0.1) up to a 5% and 10%, respectively. 
For both cases the fits found in the SC model approach are 
always less than 1% away from to the exact behavior. 

Fig U displays the ^-dependence (combining both the 
strong and weak dependences) in the skewness and kurtosis 
of the velocity divergence for different values of the spectral 
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where B® = (5q) /(<5q) corresponds to the amplitude of 
the (dimensional) skewness of the initial density field. 
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Figure 8. Values for the skewness and kurtosis of the velocity 
divergence at tree-level for open models. The smoothing effects 
strongly suppress the Q-dependence, which increases with the or- 
der of the hierarchical amplitudes Tj. The dotted line displays the 
n = —3 (unsmoothed) case while lower lines depict the smoothed 
cases for n = —2, —1, 0, 1, respectively, from top to bottom. 



index. There are two main features that arise in these depen- 
dences. First, although the relative variation of Tj with SI 
increases with the smoothing effects (larger spectral index), 
the absolute variation decreases with them, what makes it 
more plausible to probe S7 from lower values of the spec- 
tral index (which effectively means looking at small scales). 
Secondly, the dependence with the density parameter in- 
creases with J, so that looking at higher-order moments of 
the velocity divergence in the galaxy catalogues seems to be 
a promising test for determining S7. 

We stress that the behavior of the closed models is just 
an extrapolation of the one seen for the open models near 
S7 = 1, so that only small changes arise in the range 10 > 
> 1, as compared to the case of the open models. 



3.3 The ^-dependence for Non-Gaussian Initial 
Conditions 

For both cosmic fields (density or velocity) the Si- 
dependence with NGIC (for dimensional scaling, {8 J ) C tx 
(S 2 ) J ^ 2 ) can be recovered by combining in a straightforward 
way the results we have presented so far. 

This can be done by noticing that the NGIC results are 
given in terms of the Gaussian tree-level amplitudes (see §2.3 
for the velocity and Paper II for the density). Thus, to obtain 
the ^-dependence one just have to substitute the Gaussian 
tree-level amplitudes (i.e. ,5^ or Tf) with the corresponding 
Sl-dependent amplitudes, given above in Eq.Ji?! or Eq.fji]]. 
For example the variance a 2 of the density field for NGIC 
is: 



of + s 2 ,3 erf 



*2.3 



'13 
,21 



+ 



_6_ 

21 



H 

rr 3/91 



-l 



B 



(37) 



4 DISCUSSION AND CONCLUSIONS 

For Gaussian initial conditions (GIC), we found in Paper 
I, that the Spherical Collapse (SC) model gives very accu- 
rate predictions for the hierarchical amplitudes Sj in the 
perturbative regime, as compared to the results derived by 
Scoccimarro & Frieman (1996a, 1996b) for the exact PT 
in the diagrammatic approach (see Paper I §5.1). We also 
compared the predictions for the higher-order moments from 
the SC model to those measured in CDM and APM-like N- 
body simulations, and they turned out to be in very good 
agreement in all cases up to the scales where ai « 1, sup- 
porting our view that the tidal effects only have a marginal 
contribution to the reduced cumulants (see Paper I §5.2 & 
§5.3). For Non-Gaussian initial conditions (NGIC) the SC 
recovers all terms given by the exact PT for the variance, 
a, skewness, S3, and kurtosis, S4, up to non-local integrals, 
involving initial n-point functions that arise as a result of 
the coupling between the asymmetric initial conditions with 
the tidal forces (see Paper II §3.4). The measured higher- 
order moments in the N-body simulations with NGIC (with 
dimensional, texture-like, scaling) turned out to be in good 
agreement with predictions from the SC model (see Paper 
II §3.3), what gives further support to the domain of appli- 
cability of the SC model in PT. 

In this paper, we have extended the work presented in 
Papers I & II and have applied the SC model in Lagrangian 
space to study the predictions for the cumulants of the ve- 
locity divergence field, as well as the Sl-dependence of the 
one-point cumulants in PT. We provide with the results for 
the velocity field up to one loop (first correction to the lead- 
ing term). However, the loop predictions turned out not to 
be as accurate as those for the density field, suggesting they 
are dominated by tidal effects. This can be explicitly seen 
in the variance, see (Eqs.[^|-[Q) where the first non-linear 
correction to the linear contribution is only accurate within 
an order of magnitude with respect to the exact result for 
the unsmoothed field. The situation might be different for 
the smoothed cumulants, as we have seen (in Papers I&II) 
that for n ~ —1.5 tidal effects seem to cancel out, at least 
in the case of the density. The SC loop predictions for the 
reduced cumulants, are in better agreement with those ob- 
tained in the exact PT as comparison with the results of 
Scocimarro & Frimean (1996a) and Protogeros & Scherrer 
(1997) show. The latter suggests a cancellation of the tidal 
contributions in line with the behavior found for the density 
field. 

As for the S7-dependence, we showed that the SC model 
is the exact solution at tree-level for the one-point cumu- 
lants for arbitrary density parameter, Q (see also §|a| for a 
demonstration in the context of B92). This property was 
explicitly used to derive the S7-dependence of the cumulants 
at tree-level for the density and velocity divergence fields. 
As an independent check to our results we worked out the 
spherically symmetric solution to Lagrangian PT, which was 
shown to be equivalent to our solution within the SC model 
in Lagrangian space as well (see §[b]). This way we were able 
to compare available fits to the S7-dependence of the hier- 
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archical amplitudes at tree-level in the Lagrangian PT ap- 
proach to our fits within the SC model approach. We found 
agreement within a few percent between these results as well 
as with the exact numerically integrated solution (see Figs. 
U and 7 respectively) . 

We also present, as a new result, the f2-dependence of 
the one-loop corrections in the SC approximation, for GIC 
and NGIC. We show that both the variance and the Sj 
ratios increase as decreases (fl < 1). For the variance 
with GIC, this trend is in line with fitting formulas and 
numerical simulations (Peacock & Dodds 1996) and reflects 
the slowed-down growth of perturbations on large scales for 
high density FRW universes. Although this effect is small 
(less than 10% for a < 1) the SC results should be accurate 
for n ~ —1.5, and could therefore be useful when compared 
to observations of the large scale structure in the weakly 
non- linear regime. 
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APPENDIX A: GAUSSIAN TREE-LEVEL WITH 
ARBITRARY O. 

We want to show here, in the framework of Bernardeau 
(1992) calculations, that the equations of motion that gov- 
ern the tree-level amplitudes for Gaussian initial conditions 
in a FRW universe with arbitrary density parameter, Q, are 
those of the SC dynamics. 

Let us first introduce the vertex generating function 
which is defined as (see B92), 



Gf(t) 



E 

k=l 



(-r) k 



where Vk = (-F )™ , for a generic field F and the superscript 
n in the average denotes that all tree-level contributions 
come through diagrams that have n external lines (see B92 
for details). According to B92, the equations of motion for 
Gs have the form, 



a d-a +f{a)T dT- 
d 



Gs(a,r) 
d 



a 7T +f( a ) T W 
. da or 



[1 + Gs(a,r)] GA<s>(a, t) 
1 



<?A$(a,T 



■Ga^{o,,t) = 



2v 



1 d log v 
- + 



d log a 



<5a* (a,r) — -QQs(a,r) (Al) 



being, $ the peculiar velocity potential, f(a) = 
d log D(a)/d log a, and v = dloga/dlogt = Ht, where 
t is the comoving time. First, as we want to recover the 
equations of motion in terms of the density contrast and 
the velocity divergence, we shall assume at the tree-level, 
Gs(a,r) — 8, as pointed out by B94 although there it was 
restricted to flat-space. Replacing this last identity in the 
continuity equation [0 leads to <5a* = 6 = V • V/f(a) H 
provided we introduce the comoving time derivative in the 
way, 



dt 



Substituting this into ( |Al[ ), we get 

^ +2H ^tjTT8) = ^ nH2s ( 1+S ) = ^GpS (1+6), (A3) 

where in the first equality we have made use of the equation 
of motion, H = — (1 + 0/2) H 2 valid for a A = universe, 
while the second equality follows from the Friedmann equa- 
tion. Equation |jA3j is exactly the equation that describes 
the evolution of density perturbations in the SC model (see 
Paper I §4.1). 
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On the other hand, since in the SC model in the pertur- 
bative regime the evolved density field is given by a local- 
density transformation of the kind, 



(A4) 



from the above equivalence between the vertex generating 
function and the density field from the SC model at tree- 
level, 



5{5i)=g s {a,-T), 



(A5) 



it follows that, Cfc(fi) = i>k(Sl), i.e., the coefficients of the 
local-density transformation are those of the SC model for 
a FRW universe with arbitrary Q. 



APPENDIX B: THE ^-DEPENDENCE IN 
LAGRANGIAN PERTURBATION THEORY 

In this section we give a brief account of the Lagrangian 
PT which is based on an expansion of particle trajectories 
around the initial positions. We proceed and show the for- 
mal equivalence of the latter to the SC model in Lagrangian 
space. This result along with the available fits to the La- 
grangian PT found by Bouchet et al., (1995, BCHJ95 here- 
after) allows to extend their results for the skewness, to the 
kurtosis of the smoothed density and velocity fields. 

Lagrangian PT contains the ZA as the first order of 
the perturbative series. Higher-orders were considered by 
Moutarde et al., (1991), and were later generalized by 
Bouchet et al., (1992). Our analysis here essentially follows 
the steps of BCHJ95 were the perturbative analysis was ex- 
tended up to third order. 

We start by introducing the density contrast and relat- 
ing it to the displacement field, i.e., the perturbation around 
the original trajectory in Lagrangian space, 

7-'. 



1 = 



where J is the Jacobian of the transformation between real 
(Euler) and Lagrangian space, x(t,q) = q + ty(t,q), 

dxj 
dx 



Jin = TT— = 5ij + 



J 



dq 



being ^ the component of the displacement field along the 
i direction. 

Expanding J and, consequently, ^ in a perturbative 
series, 



J = l + eJ^+e 2 J^+- 



(Bl) 



and rewriting J in terms of the perturbative contributions 
from the displacement field, we have, 

J = l+eK w +e 2 (K {2) +L {2) )+e 3 (A- (3) +L (3) +A/ (3) )+- . .(B2) 

(see also Eq.[17] in Bernardeau 1994c) where the invariant 
scalars K, L, M are, 



K = V ■ # = ^2 *i 



M = detf* 



t,ji (B3) 
form which the first perturbative orders are read off, 



K 



(m) 



"(2) 



L C3) = 



i 

- V(* (1) * (1) 

O / y y J.J l ,J ?>* 



det[*«]. 



M 



(3) 



1,3 J,» ' 



The fluid equations of motion in Lagrangian space for a 
collisionless fluid can be written down in the simple form 
(see BCHJ95 for details), 

J(T,q)Vx = 0(T)[J(r,q)-l], (B4) 

where dr oc a" 2 dt, and the dot denotes time derivatives 
with respect to r, while, 

6 



/9(r) = 



(B5) 



r 2 + k(Q)' 

with k(fl) is a function of the geometry of the universe, 
k{Q = 1) = 

k(n < i) = -l 

k{Q > 1) = +1. 

The conformal time parameter r = t -1 / 3 for a flat universe, 
and is defined as r =1 1 — | _1//2 , for non-flat universes. 



Equation [B4 has a unique solution provided the fluid is 
irrotational V x x — 0, what we shall assume in the following. 

Introducing the perturbative series in the equations of 
motion [ B4 1 , one finds as generic solutions to the perturbed 
equations of motion for an irrotational fluid, 



K 



{1) (r,q)= gi (r)K W (n,q) 



K <2> (r,q)=g 2 (r)K {2 >(r l ,q) 



(2), 



K (3) (r, q) = 53a (r) M w (n, q) + g 3b (r) L w (n, q) 



(3), 



(B6) 



The first two orders are separable, unlike the third one, 
where two independent growing modes have to be taken into 
account. The above expression shall be taken as a definition 
of the gk functions which encode the f2-dependence of the 
Jacobian J, and thus that of the equations of motion in 
Lagrangian PT. 

We now impose the spherical symmetry in the growth 
of perturbations which, in terms of the displacement field 
reads, 



=0 for i + j. 



(B7) 



Applying these simple properties on the perturbative series 
of the Jacobian we find, 



A- (1) (r,q) 
K i2) {T,q) 
L i2) {r,q) 



3*«(r,g)=3 5l (r)*^(r l ,g) 



Sgi(r)^l(n,q) 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological Perturbation Theory and the Spherical Collapse model 13 



K (3) (r, q ) = (gza(r) + gs b (T))^ (n,q) 



L^(r,q) = 6 52 (T) Sa (T)^y(n,<z) 



Ci)' 



M (3) (r,q) = 9?(r) 



(i)' 



(B8) 



where we have made use of the fact that K^ 2 \Ti,q) = 
L {2) (Ti,q) (see BCHJ95) which yields ^ k {n,q) = 

^kk ( T i^l)- Replacing these last expressions in the pertur- 
bative orders of the Jacobian one gets, 



J w (r,g) = 3pi(r)*W(r 4 ,«) 

J (2) (r,q) = 3[g 2 1 (r)+g2(r)]¥ k 1 l 2 (T 1 ,q) 

J (3) (r,q) = [g 3a (r)+6g 3b (r) 

+ 6g2(r)gi(r)+g 3 1 (r))]^[ 1 l 3 (n,q). (B9) 

On the other hand, the unsmoothed transformation of 
the evolved density contrast in the SC model reads, 

X X i V2 X 2 i ^ 3 X$ i 

o = (>i + y oi + — Si H 

while the analog expansion in terms of the Jacobian gives, 
5 = ~eJ W +e 2 (J {1)2 - J (2) ) 

+ e 3 ( _ j( l)« + 2J (l) J( 2)_ J( 3) )+ ... 

from which we see that the linear term is just <5; = —e 
and all perturbative orders can be expressed as powers of the 
linear term, since J (n) oc [J (1) ]", in the same way <5 (n) oc [5i] n 
in the SC model. 

Identifying in the same way the higher-orders one by 
one (in powers of e) and replacing the expressions for the 
J (i) according to Eq.pfl]. we arrive at, 



1>:\ 



10 _l2^ + ^ + 6^ 
91 91 91 



(BIO) 



Making use of these coefficients we immediately obtain the 
tree-level amplitudes for the skewness and kurtosis of the 
density field in terms of the g k functions which, in turn, 
incorporate the ^-dependence, 



1 
9 



92 



272 + 150 7 + 21 7" -12 (24 + 77) „ 

91 



+ 48 I ^ +8^+48^ 
Ji J 91 91 



(B12) 



where the ZA (or the SSZA in our spherically symmetric 
approach), which corresponds to the first perturbative order 
in the Lagrangian approach, is recovered by setting gk = 0, 
i.e., the ^-independent term. 

Though the exact analytic expressions for (72/91 are 
known (see BCHJ95), no exact results are available for 
93a /gf and g3b/gi, so we shall focus on the numerical fits for 
the cosmologically favored parameter space for f2 and work 
out the predictions. According to BCHJ95, near 0,-1 good 
fits are provided by the factors given by Eq.j3l]]. Replacing 
these expressions for the ^-dependence in the skewness and 
kurtosis, we finally get the fits to the smoothed Sj ampli- 
tudes given by Eq.[B2j. 

Similarly, we can derive the corresponding fits to the 
weak ^-dependence in Lagrangian PT for the skewness and 
kurtosis of the velocity divergence just by recalling the re- 
lationship between the unsmoothed coefficients of the local 
transformations for the density and velocity field provided 
by Eq.|Q, and the expressions for the smoothed coefficients 
in terms of the unsmoothed ones provided by Eq. Q . 

The smoothed T3 and T4 then read, 



T s (fi)/(n) 



T 4 (n)f(n) 2 = ~ 



24 + 26 7 + 7 7 2 -4 (22 + 137) ^ 

9i 



+ 64 I ^) +8^+48^ 
Jl 91 91 



(B13) 



where the overall factor /(fi) gives the strong dependence on 
the density parameter given by Eq.j3^]. Introducing the fits 
given by BCHJ95 for the density field in the above given 
expressions, we obtain their approximated ^-dependences 
given by Eq.pl]. 



53 = 3 ^2 = 4- 2^ 

91 

5 4 = 4i/ 3 + 12i/f 

= § L-36^+6 fcV + ^ + 6^V(Bll) 

9 ^ g\ \g\) g{ g\ J 

The expression for S3 was already derived by BCHJ95 so 
the above expressions extend their results to the kurtosis in 
a simple fashion. An immediate extension of these results to 
the smoothed fields is achieved by applying the expressions 
for the smoothed coefficients of the local non-linear trans- 
formation that describes the spherical collapse Vk (see §4.4 
in Paper I), 

S3 = 37^ = 4 + 7-2^| 
91 
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